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In the holographic model of QCD, 6 dependence sharply changes at the point of confinement- 
deconfinement phase transition. In large N QCD such a change in 8 behavior can be related to the 
breakdown of the instanton expansion at some critical temperature T c . Associating this temperature 
with confinement-deconfinement phase transition leads to the description of the latter in terms of 
dissociation of instantons into the fractionally charged instanton quarks. To elucidate this picture, 
we introduce the nonvanishing chiral condensate in the deconfining phase and assume a specific 
lagrangian for the rj' field in the confining phase. In the resulting picture the high temperature 
phase of the theory consists of the dilute gas of instantons, while the low temperature phase is 
described in terms of freely moving fractional instanton quarks. 
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INTRODUCTION 



Color confinement, spontaneous breaking of chiral symmetry, the U(l) problem and the 9 dependence are some of 
the most interesting questions in QCD. At the end of the 1970's A. M. Polyakov |l| demonstrated charge confinement 
in QEDj,. This was the first example where nontrivial dynamics was shown to be a key ingredient for confinement: 
Oh| Instantons (monopoles in 3d) play a crucial role in the dynamics of confinement in QED3. Instantons in four 
dimensional QCD were discovered more than 30 years ago @. However, their role in QCD4 is still not well understood 
due to the divergence of the instanton density for large size instantons in the confining phase. Soon after, 't Hooft 
and Mandelstam [3[ suggested a qualitative picture of how confinement could occur in QCD4. The key point in the 
't Hooft - Mandelstam approach is the assumption that dynamical monopoles exist and Bose condense. Many papers 
have been written on this subject since the original formulation Q; however, the main questions, such as, "What are 
these monopoles?"; "How do they appear in the gauge theories without Higgs fields?"; "How do they interact?"; " 

' [ What is the relation (if any) between the 't Hooft - Mandelstam monopoles and instantons?" , are still not understood 

(for a review see Q). 

We re-consider these issues from a slightly different angle by analyzing phase transitions as a function of temperature, 
T, at nonzero 9 parameter. We study the evolution of the most important field configurations as the phase transition 
' line is crossed. Indeed, understanding 9 dependence gives a very good idea about the dynamics of the most important 
color fluctuations with nontrivial topology. On the other hand, we will see that 9 dependence can be studied using 
!LJ 1 an effective lagrangian approach with color singlet degrees of freedom. 

The main result of our work can be formulated as follows. In the holographic model of QCD [|| confinement 
-deconfinement phase transition happens precisely at the value of temperature T = T c where 9 dependence experiences 
a sudden change in behavior from 9/N in the low-temperature (confining) phase to cos 9 in the high temperature 
(deconfining) phase Consider now QCD with large number of colors N. For very high temperatures T » 
T c ~ Aqcd the typical size of instantons is very small and the instanton gas is dilute with density of order e~ N . 
Calculations in this region are under complete field theoretic control and the vacuum energy behaves like e~ N cos 9. As 
the temperature is lowered to be of order T c the average instanton grows in size and the perturbative expansion around 
the instanton field configuration becomes unreliable. However, we argue, based on some reasonable assumptions, that 
in the large N limit the average distance between the instantons remains much larger then the instanton size all the 
way down to the critical value of temperature T — T c . Below T c the instanton expansion breaks down and a consistent 
field theoretic calculation with non-overlapping instantons is no longer possible. It is then natural to assume that at 
T = T c there is a sharp transition in 9 behavior, which can be associated with confinement-deconfinement transition, 
just as in the holographic model. The value of T c can be estimated by a one-instanton calculation. 

To elucidate the physics of the transition we consider a model where the chiral condensate does not vanish in the 
deconfining phase. The holographic model of QCD is a good example where this phenomenon occurs. On the field 
theoretic side this can be achieved by coupling fundamental matter to the hidden gauge group whose dynamically 
generated energy scale is higher than that of QCD. In the presence of nonvanishing chiral condensate, its phase <p 
(which can be canonically normalized to yield rf) is a perfect probe of the topological charges of the constituents on 
both sides of the phase transition line. This is a consequence of uniqueness of rf meson: it always enters the system 
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in combination (6 — ip) irrespectively whether it is in the confining or the deconfining phase. One should remark here 
that our results do not really depend on the value of Nf as long as Nf <C N. 

The plan of the paper is as follows. We start in Section II by reviewing recent work on the holographic model 
of QCD where we note that 9 behavior sharply changes at the point of confinement-deconfinement phase transition. 
We return from holographic model to quantum field theory in section III where we argue that instanton expansion 
breaks down sharply at some critical temperature T c and estimate its value in terms of Aqcd- Sections IV and V 
are devoted to the physical interpretation of the phase transition. Here we attempt to answer the following question: 
what happens to the well-defined objects (instantons) as the phase transition line at T = T c is crossed from above. 
We assume a certain lagrangian for the low temperature phase and show, that under this assumption, instantons 
do not completely disappear from the system but rather dissociate into the instanton quarks 1 , the objects with 
fractional topological charges ±l/iV which become the dominant quasi-particles. At non-zero temperatures T < T c , 
the instanton quarks carry, along with fractional topological charges, the fractional magnetic charges which makes 
them the perfect candidates to serve as the dynamical magnetic monopoles, the crucial element of the standard 't 
Hooft and Mandelstam picture for the confinement 0], see the end of section V for details. In section VI we formulate 
our main results and discuss future directions. 



II. HOLOGRAPHIC MODEL OF QCD 

In this section we consider physics of holographic model of large N QCD with Nf -C N flavors [f| @ . This section 
is the review of previous work; our main point here is the observation that 9 dependence changes once we go from the 
confining to the deconfining phase. The holographic model of QCD is realized by placing Nf D8 — D8 pairs in the 
background created by N DA branes. In the weak coupling regime the D4 branes span x° . . . x 4 coordinates, while 
D8 — D8 branes are pointlike in the x 4 direction and span the rest of spacetime. One of the directions (denoted by 
X4 below) along the DA branes is compactified on a circle of radius R4 (which sets the scale of the glueball masses) 
with antiperiodic boundary conditions for fermions. The value of the asymptotic separation between the D8 and D8 
branes, denoted by L, is a parameter of the brane construction along with R±,N,Nf and string coupling and length 
g s and Z s (which will be set to unity in the rest of the paper). It will be convenient to introduce the five-dimensional 
t'Hooft coupling A5 = g s Nl s and its four-dimensional counterpart defined at Kaluza-Klein scale A4 = A5/-R4. In 
the limit A4 -C 1, Aqcd *C I/R4 and hence the theory approximates QCD pretty well. String theory is solvable in 
the opposite regime, A4 ^> 1, where there is no clear separation between the QCD scale and the supergravity/DBI 
dynamics. It is this regime that we consider below, in the hope of drawing some qualitative lessons. 

For A 4 3> 1 the rules of gauge/string duality [l(| HH, 03 instruct us to pass to the metric which is the product of 
the D4 branes backreaction: 

ds 2 = (|) 5 ({dx,f + f(u)(d X r) + (I) * (7^ + u2 ^i) (i) 

where f(U) = 1 — U^/lf 3 and the U coordinate is bounded from below by Uk- The (U, x 4 ) are the analogs of polar 
coordinates on the plane, which is ensured by the relation 

= (2) 
3 \Uk) 3 \Uk) y ' 

where in the second equality we used the relation between the t'Hooft coupling and the curvature scale of the space 
©■ 

As explained in [131 ] , the inclusion of the 9 angle in this model corresponds to having nonvanishing integral of RR 
one-form over the x A circle 

/ Ci = / F 2 = 9 mod 2tt/c; F 2 = dd (3) 
Js 1 Jd 



Instanton quarks originally appeared in 2d models. Namely, using an exact accounting and resummation of the n-instanton solutions 
in 2d CP N ~^ models, the original problem of a statistical instanton ensemble was mapped unto a 2d-Coulomb Gas (CG) system of 
pseudo-particles with fractional topological charges ~ 1/iVjSj. This picture leads to the elegant explanation of the confinement phase 
and other important properties of the 2d CP N ~ 1 models [8j. Unfortunately, similar calculations in id gauge theories is proven to be 
much more difficult to carry out 
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where k is an integer number and the integral is over the S 1 parameterized by x 4 e [0, 27ri?4]. In the first equality 
we used the fact that the (U, x 4 ) space has the disk topology and Stokes theorem. One can solve the equation of 
motion for F2 without taking back-reaction into account (which is justified as long as Nf <C N) and substitute into 
the action; the result for the vacuum energy at small 9 is [13| 

E vac ^^-e 2 (4) 

where \g ~ i s the topological succeptibility. The addition of fundamental matter results @ in the effective 

lagrangian consistent with Veneziano-Witten formula for the rf mass: 2 

£e// = r(3 M 7/) 2 + — x 9 - + --M , (5) 



2 V p " 2 AS \N N f v 

where we included some numerical factors such as y/2 and yj Nf into the definition of f v i to simplify notations in the 
following sections. This result is not significantly changed when the finite temperature is introduced, as long as the 
theory is in the confining phase and the topology of the space remains the same. Eq. (jU) is consistent with the fact 
that in the confining phase physics is expected to depend on 9 via the combination 0/JV[l3j|. 

E vac = N 2 mmh^^-y (6) 

where h(x) is some function which satisfies h(Q) — h'(Q) — 0. Eq.® can also be understood from QFT viewpoint 
for finite TV as a result of summation over different branches in pure SU(N) gluodynamics, see section III of ref.[l8j 
where connection with approach [l3| is discussed. 

As we will see below, instantons are not well-defined objects in this phase. Indeed, this would contradict 6/N 
dependence since each instanton comes with an integer multiple of 9. In the holographic model this is resolved by 
identifying instantons with euclidean DO branes wrapping around the x direction which tend to shrink to zero size 
and disappear 

At finite temperature the model exhibits confinement/deconfinement and chiral phase transitions fl9L [2p| . Two 
possible metrics with euclidean time £g compactified on a circle with circumference [3 are ([!} and its double analytic 
continuation, 

ds2 = (£) 2 ( {dXl)2 + Mu)dtl + {dxi ^ + (I) * (jS) + U2dnl ) (?) 

where fa = 1 — U^/U 3 and /3 = 4^ (^{f^J 2 ■ Since the two metrics are the same, the comparison of the free 

energies is simple: as soon as Ut > Uk the black hole metric ([7]) becomes preferred. This corresponds to confine- 
ment/deconfinement transition at T = l/27ri?4. The Polyakov loop, which is the order parameter for confinement, 
vanishes in the confining phase ([I]) and has a non- vanishing value in the deconfining phase |7]). In the deconfining 
phase the x 4 circle does not shrink to zero size and Stokes theorem makes it possible to have vanishing F2, which 
minimizes the energy Q. That is, in this phase it is possible to have 

C x = ^-dx 4 (8) 

Z7T 

This leads to \g — to order N°; this is also consistent with the fact that instantons are well defined objects in this 
phase, and come with the factor of e m6 . In the holographic model this is again a consequence of the topology in the 
deconfincd phase, where the DO brane wrapping the x 4 circle cannot shrink to zero size and disappear. The factor of 
e ln in the DO brane action follows from (jSJ). Hence, we observe that the 9 dependence is different in the confining 
and deconfining phases. We will also see that such change in the behavior is also supported by analyzing instantons 
in field theory, see next sections. 

Another comment we would like to make is the existence of the phase where the glue is deconfined, but chiral 
symmetry is broken. While it is not necessarily true that such a phase exists in QCD (after all, the holographic model 
contains two variable parameters, as opposed to Aqcd), we discuss a field theoretic model with this property to 
illuminate the topological charges of the relevant constituents in the confining (section V) and the deconfining phases 
(section IV). 



2 Earlier works on the holographic derivation of lagrangian include and |15| |. Theta dependence in the holographic models has 
been also considered in Il6t LL3l . 
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III. LARGE N QCD AT T > T c 



In this section we estimate the value of T c where the instanton expansion breaks down and the 9 dependence 
presumably experiences a sharp change. In the regime T > T c the 9 dependence is determined by the dilute instanton 
gas approximation (See below for the discussion of the assumptions that are necessary for this statement to hold). 
Applicability of the instanton expansion (small density) implies that the 9 parameter enters partition function in a 
very simple way ~ exp (— iff). The instanton contribution to the rf mass ~ exp (—j(T)N) is exponentially suppressed 
for any small (but finite) positive 7 > 0. In contrast: at arbitrary small and negative 7 < the instanton expansion 
obviously breaks down. The 9 behavior presumably drastically changes at T c determined by 

7 (T = T c ) = ==► T c = cA QCD . (9) 



A. Instantons at T > T c with (0|^|0) ^ 



In the following we will be interested in the instanton density in the dilute gas regime at T > T c . We assume that 
the non-vanishing chiral condensate (0|t/>^|0) exists in this region. 

The instanton-induced effective action for Nf massless fermions can be easily constructed. In particular, for Nf = 2 
flavors, u, d the corresponding expression takes the following form, [2ll I22L I23I Izj] . 



3 

32 



dpn(p)[ ^ttV 



(u R X a u L )(d R X a d L 



{uRU L )(d R d L ) + 
3 

-r{uR(J^X a u L )(d R a^X a d L ) 



(10) 



H.c. 



We wish to study this problem at nonzero temperature and small chemical potential p (to be discussed later in the 
text) for T > T c , and we use the standard formula for the instanton density at two-loop order [25| 



n(p) = C N {fii{p)Y N p- b eM-Pn{p)\ x exp[-(A7^ + ^{2N + N f )^T z )p% 



(11) 



where 



Cn 



0466e^^L34^ 
(JV-l)!(JV-2)! 



b 



-b\og( P A QCD ), (} n {p) = mp) + Yb log (nr* 
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b' 



34 



N 



—NfN+ — . 
3 N 



3 3 11 3 

This formula contains, of course, the standard instanton classical action exp(— 87r 2 /<? 2 (p)) ~ exp[— (3i(p)] which 
however is hidden as it is expressed in terms of Aqcu rather than in terms of coupling constant g 2 (p) . By taking the 
average of eq. lfTTJ)) over the state with nonzero vacuum expectation value for the chiral condensate (0|^|0) 7^ 0, one 
finds the following expression for the instanton induced interaction, defined as Vinst^) = — (0|iinst(<^)|0), 



(0|V^|0) iV/ cos(^ - (9) 



dpn(p)[ ^tt 2 p 3 



2 



Nf — l 



~ a ■ A QCD c °s(<y5 - 8), 



(12) 



where we introduced small dimensionless parameter a <C 1 which essentially governs the relevant physics. We assumed 
factorization for the chiral condensates in large N limit in deriving (TT2")) . Hence the square bracket in the eq. (flU)) 
vanishes. We also assumed that the condensates for all flavors are equal, (0\u r ul\0} — (0\dRdz\0) = ... = 2(0|^f|0). 
Finally we introduced a singlet phase of the chiral condensate (f(x) which we iddentify with the rf{x) field in the 
standard way, 'tpRtpL ~ e lv ^ x ^ Nf . As expected, the 77' field ip(x) enters the lagrangian in unique combination with 9 
as (tp — 9) which is consequence of the anomalous Ward Identity 3 . 



In the presence of the massless chiral fermions the 9 dependence as is known goes away in full QCD. To avoid any confusions later in 
the text we remark here that our discussions of 8 dependence in this paper deals exclusively with the dynamics of gluons when light 
fermion degrees of freedom are frozen such that essentially we analyze the 9 dependence in gluodynamics rather than in full QCD. In 
different words, we assume a quenched approximation for Nf <C N. Precisely the 9 dependence in quenched approximation plays a 
crucial role in understanding of the dynamics of strongly interacting systems. 
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The mass of the <p field in the chiral limit is determined by the instanton density in this phase and is expressed in 
terms of parameter a, 

L = -\f^^) 2 + mfrft cos(^ - 9), m%, f* = a ■ A 4 QCDl (13) 

where f n i is defined in the standard way as a normalization of <p field, rf{x) = f^/ipix), and in general f v '(T) depends 
on temperature T (though it will not be explicitly shown later in the text). We also keep only the lowest Matsubara 
frequency for ip field in the environment with T ^ which ensures the validity of the static approximation for all 
interactions involving ip. This effective lagrangian, is by definition a Wilson type lagrangian for the light rf field which 
is valid as long as rf field is light, 

m v > ~ \/~ak QC D < -Agcn, (14) 

In the large N limit parameter a ~ e~^ N is exponentially suppressed 4 for temperatures above T c , a -C 1 and the 
instanton expansion converges. For T < T c the instanton expansion makes no sense (breaks down) and the expansion 
parameter becomes large a 3> 1. We assume that 9 dependence sharply changes at T — T c . We estimate the value 
of T c by equalizing 7 = according to eq. §§§ see below. 5 In deriving the low energy effective lagrangian for the rf 
field we should, in principle, use the exact formula for the instanton density and not pip which is only valid in the 
two-loop approximation. We assume that the perturbative corrections for T ~ T c , although large, do not drastically 
change the physics. Then we will see that for any T > T c the dilute instanton approximation is valid, since the 
average distance between the instantons is paramctrically larger then their size, sec cq. (|19[) below. To reiterate, we do 
not know how to do an honest instanton calculation in the close vicinity of T c , but we assume that the perturbative 
expansion around the instanton field configuration can still be performed and would yield a ~ e ~~<( T ) N where j(T) 
is a monotonic function vanishing at T = T c . Then, for T > T c the dilute instanton gas approximation is good, for 
T < T c it is no longer valid, while T — T c describes the phase transition point with drastic changes in 9 behavior. 

We should also note that one can estimate T c (fj) for non zero chemical potential /i ^ as long as the chiral 
condensate does not drastically varies with fi, which we assume to be the case at least for sufficiently small /i. It 
allows us to estimate not only a single point T c on the phase diagram but entire phase transition line T c (p) for 
sufficiently small /i <C T c . 



B. Numerical estimates 



First, we estimate the critical temperature T c by solving eq. ([9]) and calculating coefficient c using the expression 
for the instanton density (|11[) . As the first approximation (which greatly simplifies computations) we neglect all 
log(pAQC£>) factors in evaluating J dp integral. In this case the integral can be computed analytically and the limit 
N — > oo can be easily evaluated. The result for the instanton contribution takes the following form (as expected) 
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nT 



— In 

3 V Aqcd 



- 1. 



(15) 



where we neglected all powers N p in front of e lN (as it does not have any impact on computation of T c at N = oo) 
and used the standard Stirling formula 



(16) 



4 See also [2(1 for earlier discussions on the subject. 

5 It is conceivable that the phase transition and sudden change in 9 behavior occur at the same point T c for any finite N, and not only for 
N = oo. This assumption allows us to make some reasonable estimate for T c for finite N. By obvious reasons, an estimate of T c at finite 
N suffers from some inherent uncertainties. Indeed, T c in this case is determined by an approximate condition a ~ 1 in contrast with 
precise equation ||9} valid for TV = oo case. The condition a <JC 1 implies that the r/ field is much lighter than all other degrees of freedom 
in the system in the chiral limit and condition (1141 is satisfied. It is clear that this condition can be always satisfied for sufficiently large 
N where parameter a is exponentially small at T > T c . When T becomes close to T c from above, parameter a increases and becomes 
order of unity at some point. This is precisely the region where instanton approximation breaks down. Therefore, according to our logic, 
the 8 dependence may sharply change here. We identify this point where a ~ 1 with the point of the phase transition T c . Of course we 

do not know the precise coefficient here (magnitude of a could be, for example 3, instead of I), but the extracting of a large power in 

_ 3 

such an estimate , T c ~ P^qcd ■ a 11N should not produce a large error for estimation of T c even for physically relevant case N = 3. 
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to evaluate N — > oo limit. 

There are three main reasons for a generic structure (|15p to emerge. First of all, there is an exponentially large 
"T— independent" contribution, expressed as e +186JV in eq. (fT5| . This term basically describes the entropy of the 
configuration such as number of embedding SU(2) into SU(N) etc. Secondly, there is a "T— dependent" contribution 
to Vinst(y) which comes from J n(p)dp integration (jlip . It is proportional to 



A 



QCD 



= cxp 



AT • In 

3 UqcW 



(17) 



Finally, all fermion related contributions such as a chiral condensate or non-vanishing mass term enter the instanton 
density as follows - (O^O)^ - e JV.(«lu|(o|^|o>|) i For k = -^f — ► this term obviously leads to a sub leading 
effects 1/N in comparison with two main terms in the exponent (|15p . Therefore, such terms can be neglected as 
they do not change any estimates at A^ = oo. It is in accordance with the general arguments suggesting that the 
fundamental fermions can not change the dynamics of the relevant gluon configurations as long as Nf <C N. Indeed, 
the formula for j(T) for pure gluodynamics is given by the same expression (fTS)) as it should be for Nt <C N. If 
the chiral condensate vanishes, one can replace it by a small but nonzero value for the quark's mass to proceed with 
our calculations. It would not alter the equation (TT5)) . Therefore our estimate below (TT5)) is not effected whether the 
chiral condensate develops or not. In different words, we essentially study a pure gluodynamics. Our treatment of 
the problem is equivalent of a quenched approximation for Nf <C AT. The fermion fields in the present study play 
an auxiliary (not a dynamical) role as a probe of the topological charges relevant for the phase transition as will 
be discussed in the next sections. When a number of fermions increases and Nf ~ N we can not proceed with the 
estimations as we have done above. In this case we do anticipate a strong dependence on fermion properties such 
as quark's masses and the chiral condensate as argued in recent paper [27] . Lattice simulations also suggest that for 
physical values for the quark's masses one should expect a smooth crossover rather than a first order phase transition, 
see e.g. [28, 29j. We remind the reader that for pure gluodynamics which is the main subject of the present paper for 
all AT > 3 and small number of flavors Nf -C N one should expect the first order phase transition. 

From discussions above it should be obvious that there will be always a point T c where two leading contributions 
with exponential e N dependence cancel each other. As a result, at A~ — * oo for T > T c the instanton gas is dilute 
with density e~ lN , 7 > which ensures a nice cos# dependence, while for T < T c the expansion breaks down, and 9 
dependence must sharply change at T < T c . We identify such sharp changes with first order phase transition. 

As explained above, the critical temperature is determined by condition 7 = 0. Numerically at one loop level 
approximation, it happens at 



7 = 



3 (aqod 



= => T C (N = 00) ~ 0.53Aqcd, (18) 



where Aqcd is defined in the Pauli -Villars scheme. A few remarks are in order. 

a. Our computations are carried out in the regime where the instanton density ~ exp(— jN) is parametrically 
suppressed at A~ = 00. From eq. (|15|) one can obtain the following expression for instanton density in vicinity of 

T>T C , 

a~cos(v~9)-e~ aN ( 1 ^), 1 » (^jr^j » W- (19) 

where a is a numerical coefficient of order one. Such a behavior does imply that the dilute gas approximation is 
justified even in close vicinity of T c as long as T Z Tc -k. In this case the diluteness parameter remains small 6 even 
in the close vicinity of T c . Therefore, the 9 dependence, which is sensitive to the topological fluctuations only, remains 
unaffected all the way down to the temperatures very close to the phase transition point, T = T c + 0(1/N). We can 
not rule out, of course, the possibility that the perturbative corrections may change our numerical estimate for T c . 
However, we expect that a qualitative picture of the phase transition advocated in this paper remains unaffected as 
a result of these corrections in dilute gas regime. 

b. In our estimate for T c we neglected (log ph.QCD) k in evaluating of the J dp integral. One can easily take into account 
the corresponding contribution by notice that J dp is saturating at p ~ (7rT) -1 . The corresponding correction changes 



This should be contrasted with the the standard requirement for finite N when the condition a ~ {A-QCD /T) b >C 1 can be only achieved 
when the temperature is very large, T 2> ^QCD- For large N the condition o < 1 is satisfied as long as T ~ Tc 3> as can be seen 
from eq. 11191 . 
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our estimate (fT8|) very slightly, and it will be ignored in what follows. Numerical smallness of correction is due to the 
strong cancellation between the second loop contribution in the exponent (term proportional to b' /b) and the first 
loop contribution in the pre-exponent in eq. (fTTjl . 

c. The transition to a different scheme leads to very large changes in the instanton density. For example, transition 
to the so-called MS -scheme is achieved by replacing e _1679Ar in the expression for Cjv, see eq. pr]) . as follows 
e -i.679JV _^ g (-i. 679+3. 721)at w jth a number of other changes, see e.g. (23|. The corresponding results would be expressed 
in terms of Aq^ d , where MS stands for MS -scheme, to be distinguished from Aqcd which is defined in the Pauli 
-Villars scheme and will be used through this paper. We shall not elaborate on these numerical issues in the present 
work. 

d. Unfortunately, we can not compare our calculations with the precise lattice results [30( for the ratio T c j ^fa at 
large N as we compute T c in de-confined phase where the string tension a vanishes. 

e. As expected, the result (|18l) docs not depend on a number of flavors Nf nor does it depend on the magnitude of 
the chiral condensate in N = oo limit as our treatment of the problem corresponds essentially pure YM computations. 

f. For finite but large N 3> Nf the corresponding numerical estimates for T c can also be given. It can be estimated 
from condition a ~ 1. However, numerical estimates in this case would depend on the value of the E/a(1) condensate 
a (O|( - 0'0) Ar/ 10) which is not well-known for T > T c . Therefore, we shall not discuss the corresponding numerical 
estimates in the present work. 

g. A similar procedure for estimation of the critical chemical potential /t c for confinement -de-confinement phase 
transition at finite N, Nf at T ~ has been previously used in re f. 13 ll where the analogous arguments on drastic 
changes of 9 at fi = fx c have been presented, see also a review paper [32|. 

h. Once T c is fixed one can compute the entire line of the phase transition T c (/t) for relatively small fj, <C T c for large 
but finite N 3> Nf. Indeed, in the weak coupling regime at T > T c the /i dependence of the instanton density is 
determined by a simple insertion ~ exp[— Nffi 2 p 2 ] in the expression for the density pip . In the leading loop order 
Tc(fJ-) varies as follows, 



T c (n) = T c {p = 0) 



1 



3Nf/i 2 



4Ntt 2 T^(h = 0) 



/t < ttT c , Nf < N. 



(20) 



As expected, /x dependence goes away in large N limit in agreement with general large TV arguments [33j. This 
formula is in excellent agreement with numerical computations [34], [35], [36[ which show very little changes of the 
critical temperature T c with \x for sufficiently small chemical potential. In particular, even for the case Nf = 2, N = 3 
where the expression (|20[) is not expected to give a good numerical estimate, it still works amazingly well even for 
N = 3. Indeed, the result quoted in [34[ can be written as 



T c (nY at = T c ([i = 0y at 1 - 0.500(38) 



'tt 2 T c 2 (^ = 0)' at 

It should be compared with our theoretical prediction (|2"0|) for this case 



Nf = 2, N = 3. 



T c (p) tfl = T c (p - 0) 



Hi 



2w 2 T 2 (ii = 0) th 



i. It is naturally to expect that the phase transition line T c (/i) at /i <C T c from (|20p connects with the phase transition 
line at very large /i c ~ \f~N as estimated in a recent paper [27j . 



fi c {T) ~ Mc (T = 0) 



Nn 2 T 2 



3N f tf(T = 0) 



NT < /ic, Nf < N, 



(21) 



where /t c (T = 0) ~ 1.4 • Aqcd^J jfj at Nf <C N [23|. This expectation is motivated by the observation that the 

nature for the phase transition along the entire line is one and the same: it is drastic changes of 9 dependence when 
the phase transition line is crossed. Therefore, we believe that the entire line is the first order phase transition as long 
as N f < N. 



IV. DUAL REPRESENTATION 



The main goal of this section is to present the low energy effective lagrangian for rj field (fT5|) in the dual form. 
The rf field will play a crucial role in the following two sections. As we shall see in a moment the ry'-field is a perfect 
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probe of the glue configurations. This field will help us to investigate the topological charges of the constituents in 
both phases, and therefore it will help us to interpret the nature of the phase transition whose critical temperature 
T c was computed in the previous section. In section II we discussed a holographic model with nonvanishing chiral 
condensate. Here we consider a field theoretic model with this property. 



A. Coulomb Gas Representation: formal derivation 

The effective low energy dense-QCD Lagrangian (fT3")) is the sine-Gordon (SG) Lagrangian. Many of the special 
properties of the SG theory apply. One of these properties is the admittance of a Coulomb gas (CG) representation 
for the partition function. Although this is a four-dimensional theory at nonzero temperature T (rather than two 
dimensional, where all known exact results regarding SG model were discussed) and questions about renormalizability 
of the theory may come to mind, there are no such issues here since the effective action is a low energy one. Following 
the usual procedure for mapping a statistical CG model into the field theoretic SG model, the CG picture that arises 
from the effective low energy QCD action, Eq. (|13p . will be derived in this section. The statistical model contains 
some charges which appear due to the presence of cosine interaction in the field theory model. The physical meaning 
of these charges will be illuminated in the next section by analyzing the corresponding measure of the statistical 
ensemble. 

The mapping between the SG theory and its CG representation is well known. All we need to do is to reverse 
the derivation of SG functional representation of the CG in Ref. The partition function corresponding to the 
Lagrangian (|13[) is given by 7 

Z = J Vtpe-f d3x !S dThE = J Dip e -vrf*>I<lMv<pf e Xf<l 3 xco S (<p(x)-e)^ ^ 

where we introduced fugacity for the CG ensemble to be defined as, 

_ f Aqcd \ 



A=^^jaA$ 013 (23) 

Le is the Euclidean space Lagrangian. Leaving alone the integration over tp(x) for a moment, we expand the last 
exponent in Eq. (|22[) . represent the cosine as a sum of two exponents and perform the binomial expansion: 



g A / d 3 x cas(tp(x)—0) 




V (V2) M 

M=0 



M+\M-\ 

M±=0 



d 3 x M e^-oQ-M*.)-') . (24) 



The last sum is over all possible sets of M + positive and M_ negative charges Q a = ±1. The last line in Eq. ([24]) is 
a classical partition function of an ideal gas of M = M + + M_ identical (except for charge) particles of charges +1 
or —1 placed in an external potential given by i(9 — cp(x)). It is easy to see that (for a constant or slowly varying 
potential) the average number of these particle per unit of 3-volume (M}/Vs, i.e., the density, is equal to A. Thus 
making A small one can make the gas arbitrarily dilute, which is precisely the physical meaning of the fugacity. From 
this observation, one can immediately see that the average distance between charges Q a is A -1 / 3 . 

While 8 can indeed be viewed as an external potential for the gas (1241) . tp(x) is a dynamical variable, since it 
fluctuates as signified by the path integration in (f2"2"| . For each term in (pM)) the path integral is Gaussian and can be 
easily taken: 

■D^e-^fd'xiv^f £ ij:^ Q a{v ( Xa )-e) = e -i9^ =0 Qa e ~^^> b = oQaG(Xa - Xb)Qh . (25 ) 



7 To be precise, the path integral in Eq. i122l) should be understood as an integral over low-momentum modes of ip only. The upper limit 
of the momentum of ip is the ultraviolet cutoff of the effective Lagrangian ( 1131 . which should be taken as some scale smaller than T. 
Only tree graphs contribute to Z so there is no dependence on the precise value of the cutoff. 
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We see that, for a given configuration of charges Q a , —iip{x) is the Coulomb potential created by such distribution. 8 
The function G{x) is the solution of the three-dimensional Poisson equation with a point source (the inverse of — V 2 ): 

G(x a - x b ) = 1 . (26) 
Air\x a ~ x b \ 

Thus we obtain the dual CG representation for the partition function 



^ (A/2) A/ f . f - ie ~M -^-^> b=0 QaG(x a - Xb )Q b 

z = E wrin / d xi ■ ■ d XM e ^= oQa e v ■ ( 27 ) 

M±=0 + ' 



The two representations of the partition function (|22j) and (|27|) are equivalent. 

Note that the physical meaning of A is the the fugacity of the system with charges Q and it is proportional to a which 
is small in the regime under discussion. There are several important features of the action (|27p which should be noted. 
Firstly, the total Q charge of the configuration, Qt appears together with 9 which we kept as a free parameter, see (|27p . 
Such a dependence will play an important role in the following identification of Q charges as the topological charges, 
see below. The ^-dependence in CG representation (f2"T| gives an overall phase factor for each configuration. Finally, 
our dimensional parameters A and xm come into the expression (|27[) in the combination Xd 3 x which is nothing but the 

coefficient in front of the instanton contribution to the effective action S'inst = J dr J d 3 xVi ns t(<p) = — J d 3 xX cos(ip— 9) 
at nonzero temperature T, see eq. (|T2|) and definition of A (|23|) . 



B. Physical Interpretation 

The charges Q a were originally introduced in a rather formal manner so that the QCD effective low energy La- 
grangian can be written in the dual CG form (|27p . However, now the physical interpretation of these charges becomes 
clear: since <5 no t = J2 a Q a 1S * ne total charge and it appears in the action multiplied by 9 [see Eq. (|2"T)) ], one concludes 
that Qnet is the total topological charge of a given configuration. Indeed, in QCD the 9 parameter appears in the 
Lagrangian only in the combination with the topological charge density — i9G ^G ^ / (32n 2 ) . It is also quite obvious 
that each charge Q a in a given configuration should be identified with an integer topological charge well localized 
at the point x a . This, by definition, corresponds to a small instanton positioned at x a (to be precise, "caloron" at 
temperature which has topological charge Q = 1 and action 8ir 2 /g 2 independent of temperature, see [24], [25| 

for review) . To support this identification we note that every particle with charge Q a brings along a factor of fugacity 
A ~ a ([2"3]) which contains the classical one-instanton suppression factor exp(— 8tt 2 /g 2 (p)) in the density of instantons 
(|f 2p if one restores the instanton density in terms of coupling constant exp(— 8ir 2 /g 2 (p)) rather than directly in terms 
of Aqcd which is used in eq. (|f f [) and which is more convenient for numerical estimates. 

This identification is also supported by the following observation: every extra particle with charge Q a brings an 
additional weight e~ lB ® a to the partition function. This is certainly the most distinguishable feature of the non-zero 
topological charge configuration. 

The following hierarchy of scales exists in such an instanton ensemble for temperatures slightly higher than T c . The 
typical size of the instantons p ~ T _1 ~ Aq^ d The average distance between the instantons f = A -1 / 3 = Ag^a -1 / 3 
is much larger than both the average size of the instantons and the cutoff T _1 . The largest scale is the Debye screening 
length in the Coulomb gas, rjj — Ag^a -1 / 2 . This coincides with the static correlation length of the tp field, which 
is precisely 77' mass. It is important that the Debye screening length r£> is parametrically larger than the average 
distance between the instantons f, therefore large number of instantons can be accommodated within the volume 
determined by the Debye screening length rjj which justifies our Coulomb gas interpretation, at least in large N limit. 
In short: 

(size, p) <C (distance, r) <C (Debye, m) 

1 <g- 1 <g- 1 ( 28 ) 



8 One notices that the term a = b in the double sum (1 25 I t is dropped. This is the self-interaction of each charge. It would renormalize 
the fugacity A by a factor exp(— G(0)/(/ 2 , )). This factor should be dropped as it represents contribution of very short wavelength 
fluctuations of ip. Such fluctuations have to be cutoff at the scale 1/T. The self-energy of the charges comes from a much smaller scales 
which are already calculated and contained in a. 
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Due to this hierarchy, ensured by small n« 1, we acquire analytical control. In reality, of course, a ~ y Q £ D j ~ 

e~ N <C 1 is parametrically very small only at very large N while ( Aq £ d J < 1 could be very close to 1 from below. 

It implies that at N = 3 all scales could be numerically very close to each other. 

It is also quite interesting that, although the starting low-energy effective Lagrangian contains only a colorless field 
ip, we have ended up with a representation of the partition function in which objects carrying color (instantons, their 
interactions and distributions) can be studied. In particular, from the discussion above, one can immediately deduce 
that // and II interactions are exactly the same up to a sign and are Coulomb-like at large distances. 

This looks highly nontrivial since it has long been known that at the semiclassical level an instanton interacts 
only with anti-instantons but not with another instanton carrying a topological charge of the same sign. As we 
demonstrated above it is not true any more at the quantum level in the presense of the rj' field. Indeed, what we 
have found is that the interactions between dressed (as opposed to bare) instantons and anti-instantons after one 
takes into account their classical and quantum interactions, after integration over their all possible sizes and color 
orientations, after accounting for the interaction with the background chiral condensate must become very simple 
at large distances as explicitly described by Eq. (j27|) . It is impressive how the problem which looks so complicated 
in terms of the original bare (anti)instantons, becomes so simple in terms of the dressed (anti)instantons when all 
integrations over all possible sizes, color orientations and interactions with background fields are properly accounted 
for! 

Such a simplification of the interactions is of course due to the presence of almost massless pseudo-Goldstone boson 
rj' which couples to the topological charge. When the instanton gas becomes very dilute all semiclassical interactions 
(due to zero modes) cannot contribute much, since they fall off with distance faster then the Coulomb interaction 
mediated by rj' . On the other hand, when the instanton density increases when T is getting smaller, the Coulomb 
interaction becomes more screened and, as the Debye length becomes comparable to the inter-instanton distances, we 
lose analytical control. Based on this picture one can estimate the critical temperature T c where this transition must 
happen. It corresponds to the same condition a ~ 1 discussed previously in section III. 

We collect here the most important results of the present section based on CG representation (|27|) : 

a. Since Q nc t = ^2 a Qa is the total charge and it appears in the action multiplied be the parameter one concludes 
that Qnet is the total topological charge of a given configuration. 

b. Each charge Q a in a given configuration should be identified with an integer topological charge Q a = ±1 well 
localized at the point x a . This, by definition, corresponds to a small instanton (caloron at T ^ 0) positioned at x a . 

c. While the starting low-energy effective Lagrangian contains only a colorless field ip we have ended up with a 
representation of the partition function in which objects carrying color (the instantons) can be studied. 

d. In particular, II and II interactions (at very large distances) are exactly the same up to a sign, order g°, and 
are Coulomb-like. This is in contrast with semiclassical expressions when II interaction is zero and II interaction is 
order 1/g 2 . 

e. The very complicated picture of the bare II and II interactions becomes very simple for dressed instantons/anti- 
instantons when all integrations over all possible sizes, color orientations and interactions with background fields are 
properly accounted for. 

f. As expected, the ensemble of small p ~ T^ 1 instantons can not produce confinement because small instantons can 
not produce a correlation at arbitrary large distances which is a crucial feature of the confinement. This is in accord 
with the fact that there is no confinement in the high temperature phase. 

g. Physical interpretation of the CG representation l|27p is simple. The vj '— field being a dynamical field couples to 
the topological charge Q exactly as 6 parameter does due to the specific combination {<p(x) — 6) which appears in 
the low energy lagrangian. In the dual language the rj mass emerges as a result of Debye screening in the plasma of 
topologically charged instantons (interacting via rj Coulomb exchange ) similar to the well-known effect of generating 
the photon's mass in the ionized plasma due to the Coulomb interaction of charged particles. In our case instead 
of a conventional vector photon we are dealing with pseudo scalar rj field which receives its mass through the 
interaction with topological charges Q. Uncovering this picture (which allows us to measure the topological charges 
of constituents) was the main motivation for introducing the chiral condensate into the theory. 

h. We should also remark here that a similar picture for the instanton interactions occurs at large chemical potential 
p 3> Aqcd in deconfined, the so-called color superconducting phase [37]. In the present case T > T c the weak 
coupling regime (small instanton density ) is governed by small parameter a ~ exp(— 7JV) <C 1 while at large 
p ^> Kqcd,N = 3 case the corresponding small factor is (Aqcd/ ' <C l(37|. 
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V. SMALL T < T c : CONFINED PHASE. SPECULATIONS. 

In this section we want to speculate on the fate of the instantons when we cross the phase transition line at T = T c . 
To be more precise: we want to see if any traces of well defined instantons discussed above can be recovered. The 
instanton expansion is not justified in the strong coupling regime T < T c where the expansion coefficient becomes of 
order one, a ~ 1. 

Therefore, we do not even attempt to use instanton calculus or any other semiclassical computations in the present 
section. Instead, we present a few indirect arguments supporting the picture that the instantons do not completely 
disappear from the system when we cross the phase transition line from above, but rather dissociate into the instanton 
quarks [1, Q, the self-dual objects with fractional topological charges ±1/N which become the dominant quasi- 
particles. The arguments are not based on the semiclassical calculations, but rather on analysis of the low energy 
lagrangian written in the dual form similar to CG representation discussed in the previous section. Our proposal 
about the fate of instantons at T < T c originally derived in ref. [38| and to be reviewed for completeness in this section 
should be considered as one of the many possible outcomes. In this sense this section is very speculative in nature 
in contrast with previous sections where weak coupling regime is justified at large N and precise statements can be 
made. 

We start from the chiral Lagrangian and keep only the singlet rj field. We assume the following expression for the 
effective Lagrangian for the ip field which has a specific Sine-Gordon (SG) form, 

L v = \f v> {d^f + E vac cos (Zjp) , ftn§ = ^ (29) 

where E vac ~ iV 2 is the vacuum energy of the ground state in the chiral limit, expressed in terms of the gluon 
condensate, 

E vac = J(0|6£|0) = (0|6a s /(327r)G 2 |0) ~ iV 2 , (30) 

where we use the standard expression for the conformal anomaly of the energy -momentum tensor, 0£J. The ex- 
pression (|29|) of course satisfies the standard requirement crucial for the resolution of J7a(1) problem: the vacuum 
energy in gluodynamics depends on 9 through the combination 9/N. It also has a very specific SG structure for the 
singlet combination corresponding to the following behavior of the (2k) th derivative of the vacuum energy in pure 
gluodynamics 39], 

2fc t? ra\ r ™ „• „2 



0=0 ' ' 



<)f)2 , „-„ . \[dx i {Q{xx)...Q(.X2 h ))~{-f k , where Q = JL-G^G,*,. (31) 

i—l 

The same structure was also advocated in [40j from a different perspective. We shall not discuss any additional 
arguments supporting such Sine-Gordon structure referring to the original papers 9 . This is precisely the place where 
the term " speculation" from the title of this section, enters our analysis. One should also note that the combination 

_ E vac _ d 2 E vac (6) 

is nothing but topological susceptibility \g f° r gluodynamics in the large N limit. 

Now we want to represent the low energy lagrangian (|29p in the dual form (CG representation) to see if any traces 
from the instantons discussed at T > T c can be recognized. The effective lagrangian is obviously the color singlet 
object. Therefore, all color dynamics can not be recovered by this method. However, the topological charge is color 
singlet operator which is coupled to 8. The 9 parameter is not a dynamical field in QCD, however the ?/ field is, and 
it always enters the dynamics in combination (9 — ip) . Let us repeat again that this was the main reason to introduce 
the chiral condensate into the sysytem: it allows to study the dynamics of the topological charges. Therefore, in 
principle, the analysis of the rf field gives the information about the topological charges of the constituents. We use 
the trick (SG-CG mapping) below to attempt to answer the following question: what kind of constituents can provide 
the low energy behavior (|29l3ip ? 



One more additional argument supporting SG structure ~ cos (yjj^ 1X1 pure gluodynamics will be given later in the text. 
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We use the technique developed in the previous section and represent SG action in the dual form. Technically, it 
goes as follows: eq. ([25]) replaces expression (fT3|) discussed previously. As in (fT3j) the Sine-Gordon effective field 
theory (|29|) can be represented in terms of a classical statistical ensemble (CG representation) similar to ([27| with the 
replacements A — ► E vac , d 3 x — > d A x as we assume zero temperature T = in this phase. By repeating all previous 
steps we arrive at the following expression 

Ail 

AT=0 ' ' Q a =±i/N' 

The functional integral is trivial to perform and one arrives at the dual CG action, 



Z = J2 4^ / d* Xl ... f d^x M ■ I Dtpe-WJ^W ■ (e^^iO-W*.)-^) . ( 32 ) 

H.-f n « « I 1 / TV 7" 



Z = Yl M \M-\ / 1 ■ " " / M ' 6 E(<,=0 ' Qo=±1/w) ' e . ( 33 ) 



Af±=0 + 

where G(x a — Xb) is the 4d Green's function, 

G ( Xa - Xb) = 4«*{xl- Xb r (34) 

The fundamental difference in comparison with the previous case (|27p is that while the total charge is integer, the 
individual charges are fractional ±1/A. This is a direct consequence of the 9/N dependence in the underlying 
effective Lagrangian (f2U)) before integrating out <p fields, see eq. (I52"j) . 

A few remarks on physical interpretation of the CG representation (|33| of theory ([29]) arc in order: 

a. As before, one can identify Q nGt = J2 a Q a W1 th the total topological charge of the given configuration. 

b. Due to the 2tt periodicity of the theory, only configurations which contain an integer topological number contribute 
to the partition function. Therefore, the number of particles for each given configuration Qi with charges ~ 1/N must 
be proportional to N. 

c. Therefore, the number of integrations over d 4 Xi in CG representation exactly equals ANk, where k is integer. 
This number ANk exactly corresponds to the number of zero modes in the fc-instanton background. This is basis 
for the conjecture [38| that at low energies (large distances) the fractionally charged species, Qi — ±1/JV are the 
instanton-quarks suspected long ago [8|- 

d. For the gauge group, G the number of integrations would be equal to 4fcC2(G) where ^(G) is the quadratic 
Casimir of the gauge group (0 dependence in physical observables comes in the combination q^q^ )- This number 
4fcC2(G) exactly corresponds to the number of zero modes in the /c-instanton background for gauge group G. 

e. We do not use the weak coupling regime or instanton calculus anywhere in our arguments. Still, we recover the 
moduli space which we identify with strongly interacting instantons in the confinement phase of the theory. 

f. Role of the fugacity for this statistical ensemble plays E vac ~ N 2 . Therefore, an average distance between 
constituents is of order f ~ E V aJ 4 ~ Aq^A^ 1 / 2 which suggests that the system is very dense. It obviously implies 
that the instanton expansion makes no sense in this regime as all terms are equally important, which is in huge 
contrast with hierarchy from the previous case at T > T c , (|28l) . 



g. The Debye screening length rjj ~ m^, 1 ~ ^qcd^N ^> r is large. It means that the number of constituents 
participating in the screening is order of (r^/f) 4 ~ A 4 . 

h. According to eq. (|29| the number of instanton quarks in the spacetime box of size Aqcd should be A 2 as an 
average distance between constituents is f ~ A -1 / 2 . Each instanton contains N instanton quarks, hence the density 
of instantons should be of order NAq CD . 10 It is consistent with observation from holography, section II that any 
finite number of instantons will disappear from the system. 

A. The relation to other studies 

As we mentioned above our arguments in the present section look extremely speculative as they are not based on 
instanton calculus or any other dynamical calculations which include color degrees of freedom. Still, by analyzing 9 
dependence in deconfining phase we infer (indirectly) that some fractionally charged degrees of freedom emerge at 



In 13811 it was conjectured that these constituents (instanton quarks) are the driving force for the confinement. 
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low T < T c . We presented arguments suggesting that the corresponding fractionally charged constituents are the 
colored instanton quarks. We should remark here that the fractionally charged constituents have been discussed in 
the literature in a number of papers previously. In particular, there seems to be a close relation between instanton 
quarks and the "periodic instantons" [H, S, H], center vortices and nexuses with fractional fluxes 1/N, see e.g.(45j 
and references therein. We shall not discuss the corresponding connections in details in the present paper by referring 
to the original literature and the recent review by one of the author [32| where some comments on the correspon ding 



connections have been made. In the present work we want to make a few comments on two recent papers [46j . 47 1 
where the picture similar to the one presented in this work is advocated. 

We start with 46]. In that work the authors consider a specifically deformed SU(N) gluodynamics at T ^ 0. It has 
been shown that such a deformation supports a reliable analysis in the weak coupling regime in the confining phase. 
The results of the corresponding calculations imply that the relevant degrees of freedom in the confined phase are the 
self dual magnetic monopoles with action 4^ and topological charges Q — ±1/N which are precisely the features of 
the instanton quarks discussed above. 

In contrast, the starting point of ref.[47| is semiclassical calculations in the background of calorons[42} where the 
weak coupling regime can not be guaranteed. While the calculations are semiclassical in nature, and therefore, can 
not be trusted in the strong coupling regime, still, the corresponding analysis shows how well localized instantons with 
integer topological charges at T > T c may dissociate into the fractional constituents at T < T c , and become the key 
players in the confining phase. This is precisely the picture we are advocating in the present work based on analysis 
of sharp 9 changes at N — oo. It is impressive how complicated semiclassical calculations carried out in [47j lead 
to the expression for the vacuum energy E vac cos(^) advocated in [|(| using completely different technique. 11 Our 
technique does not allow us to make any dynamical calculations in this phase as all color degrees of freedom have been 
integrated out in the course of obtaining (129[) . In other words, we can not study the dynamics of fractionally charged 
constituents in contrast with papers [4y, [47[ 12 • However, the fact that the constituents carry fractional topological 
charge 1/N can be recovered in our approach because the color- singlet vf field enters the effective lagrangian as 
cos(-j^) and serves as a perfect probe of the topological charges of the constituents. One should also emphasize that 
the procedure of the recovering of the fractional topological charge 1 /N (which has been used here) is not based on 
the weak coupling expansion. 

Our final comment in this subsection is as follows. The main ingredient in holographic picture discussed in section II 
was DO brane wrapping around x which behaves differently in confined and deconfined phases, and correspondingly 
leads to a different behavior on opposite sides of phase transition line. Similar picture was also observed in ref . [481] 
where the authors studied the D2 brane in confined and deconfined phases to arrive to the same conclusion on sharp 
changes in 6 behavior. The topological objects (sensitive to 0) were identified as magnetic strings in ref.[48j]. These 
objects apparently have been observed in the lattice simulations [49|. 



VI. CONCLUSION 



We explore the consequences of the assumption that in the large N QCD and QCD-like theories confinement- 
deconfinement phase transition takes place at the temperature where the dilute instanton calculation breaks down, 
and dependence drastically changes. This assumption is supported by holographic and field theoretic arguments. At 
very high temperatures, T Aqcd instantons are well localized configurations with a typical size p ~ T^ 1 <C ^qqd- 
As the temperature is getting lower the instanton size becomes of order p ~ Aqcd however provided the perturbative 
corrections in the instanton background do not significantly change the picture, the instanton density remains dilute. 
Instanton expansion breaks down at T c , and for T < T c , the strong interacting regime and confinement are realized. 
Instantons are no longer well localized configurations for T < T c , but rather, in the picture of Section V, they are 
represented by N instanton quarks which can propagate far away from each other. The presence of the light field rj 
in the model is important for this picture and a specific lagrangian is assumed in Section V. The mass of rj in both 
phases in the dual picture can be thought as the Debye screening mass generated by the Coulomb interaction of the 
topological charges. 



11 Such a SG structure was a crucial element for recovering the fractional topological charges Q = ±1/7V in the confining phase using rj' 
as a probe, see section V and original discussions in |38|1 . 

12 In particular, we do not see a beautiful picture of a multi-component color Coulomb plasma with nearest-neighbor interactions in the 
Dynkin space advocated in [46l. |47H . Still, we do see the color- singlet Coulomb interaction of the fractionally charged ±1/JV constituents 
due to rj at very large distances where color already confined. 
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We have made a number of assumptions in the field theoretic analysis to arrive at the conclusion that the 6 
dependence changes sharply at some value of T c . We also used the holographic model to argue that this transition 
coincides with confinement-deconfinement phase transition. This conclusion is supported by the lattice simulations 

no, he nn m, m . 

The value of the critical temperature as a function of (sufficiently small) chemical potential TJu) is estimated in 
Section III. The obtained expression is in excellent agreement with numerical computations [U OS [36| • Finally, we 
presented the arguments that this line of the phase transition (which is the first order for large N and Nf <§; N) 
continuously transforms into the line fi c (T) studied at large /i in ref . [27| . The argument is based on the observation 
that the physical nature of the phase transition along entire line is the same: it is the drastic changes in 8 behavior 
when the phase transition line is crossed. 

It would be very interesting to see if Coulomb law between instantons can be understood holographically in the 
deconfining phase. As we mentioned in the text it is quite nontrivial that at large distances in the presence of rj' 
field the interaction between instantons and anti instantons is the same (up to sign) as the interaction between two 
instantons. 

Finally, we would like to make a short comment on relevance of the present analysis to real QCD with Nf = N = 3. 
The main subject of the present paper is pure gluodynamics at large N, and therefore one can not immediately apply 
the results of the present analysis to the real QCD with Nf = N = 3. However, we do expect that our results can be 
and should be compared with the lattice simulations for pur e gluodynamics for N > 3 and for QCD with Nf <C N 
when the first order phase transition is expected, see [30j, [50|, l5l|, l52|, |53j. For the case Nf ~ N our technique is not 
applicable as we explained in section IIIB. In this case one should expect that the properties of the phase transition 
is very sensitive to the details of the fermion matter fields. It a subject of a separate analysis which will not be 
discussed here. However, we expect that the picture of the phase transition as a transition between plasma phase 
(when the instanton quarks are in plasma state at T < T c ) and molecular phase (when the instanton quarks form a 
small instanton at T > T c ) qualitatively describes real QCD with N = 3 as the confinement in non- abelian gauge 
theories is determined by the dynamics of gluon (not quark) degrees of freedom. 
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